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Motivation 


□  computational  efficiency 

□  coarser  grids 

□  larger  time  steps 

spatial  concerns: 

□  how  well  are  gradients  captured?  (resolution  requirement) 
spatial/temporal  concerns: 

□  dispersion  and  dissipation  error 

□  better  characteristics  for  a  broader  range  of  wavenumbers 

scheme  stabilization  concerns: 

□  balancing  accuracy  with  stability 

□  artificial  dissipation,  filtering 
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Spectral  Representation 


Fourier  basis: 
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modified  wavenumber 
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Modified  Wavenumber 


errork  = 


20  10  5  3  (PPW) 


0 


goal:  proper  representation  of  derivative 

4 

Distribution  Statement  A:  Approved  for  public  release;  distribution  is  unlimited. 


Gradient  Capture  vs.  Resolution:  Single  Mode 


Solution/Derivative: 


X 


FFT: 


f  (x)  =  sin(x)  with  x E[0,  2jc  ] 


—  =  cos(x) 
ck 


Convergence: 


_ CD02 

_ CD04 

_ CD06 

--  CD04  (Tam,7pt) 
_ CDIO 

--  CD04  (B&B,  I  Ipt) 


PPW 
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Gradient  Capture  vs.  Resolution: 
Multiple  Modes 


Solution/Derivative: 


FFT: 


4  _  CD02 

f(x)=2sm(mx)  with xE[0,2jt]  CD04 

m=\ 


df_ 

dx 


4 

^  mcos(mx) 

m=l 


Convergence: 


_ CD06 

--  CD04  (Tam,7pt) 
_ CDIO 

--  CD04  (B&B,  I  Ipt) 
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Isentropic  Vortex: 
no  Stabilization 


CFLu  W~  0.01  a  =  M)  =  l 

CD04  f  1 1  Pt  B&B) 


need  to  stabilize 
accumulation  of  high 
frequency  error 
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Stabilization: 

Artificial  Dissipation  and  Filtering 


Artificial  Dissipation: 


dQ 

dt 


dE 

—  + 
dX 


2(-ir(Axr 


d2mQ 

dx1171 


damping  strongly  dependent  on  base  scheme 
(couples  with  temporal  scheme) 


Filtering: 

dQ  d  E 

dt  dX 


M 

i+D-ir 


m 


2m  ^2 m 


dx2 


m 


Q  = 


.V 


i+K-ir 


adds  consistent  amount  of  damping  to  base  scheme 
(decoupled  from  temporal  integration) 
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Density  [kg/m  ] 


Isentropic  Vortex: 
Traditional  Stabilization 


EFlO(l.O) 


1  + 


/A  v  10 

Ax \ 

Ci , 


distance  traveled:  SOOtakkiihs 


effectiveness  of  stabilization  strategy  dependent 
on  spectral  content 


CFLu[D  ~  0.01  a  =  1,4)  = 


CD04  Mlpt  B&B) 
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Damping  Characteristics: 
Growth  Factor 


D) 


-  need  scale-discriminant,  tunable  formulations 

-  strong  preservation  of  resolvable  modes 


0.4  0.6 

(k  A  x)/jt 
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Scale-Discriminant  Stabilization 


20  10  5  3  (PPW) 

goal:  io° 

preserve  accurately  resolved 
frequencies  and  remove  error-  10- 
prone  content 

scale-discriminant  dissipation  < —  10 

k-  k  ,  lo- 

error^  = - =  0.01 


for  the  isentropic  vortex... 

□  maintain  coherence  of  vortex 

□  preserve  vortex  intensity 


FD  stencil 

/jt 

-PPW 
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CD06 
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CD10 
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5 

CD04  (7pt  Tam ) 

0.48 

5 

CD04  (1  lpt  B&B) 

0.55 

4 
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(kA  y)/jt  (kA  y)/jt 


Isentropic  Vortex: 

Scale- Discriminant  Stabilization 


CFLuW~  0.01  a  =  M>  =  l 

CD04  f  llpt  B&B)  +  IF10  (2/3 Pi) 


ref 

IF  10  (2/3 Pi) 

EFIO 


1 6  points  across  vortex 


1.01 


0  5  10  15  20  25 

X 

distance  traveled:  SOOMAddlibs 

-  traditional  dissipation  is  sufficient  when  spectrum  is  well 
resolved 
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Isentropic  Vortex: 

Scale- Discriminant  Stabilization 


CFLuW~  0.01  a  =  M>  =  l 

CD04  f  llpt  B&B)  +  IF10  (2/3 Pi) 


ref 

IF  10  (2/3 Pi) 

EFIO 


(kA  x)/jt 


i. 

-0.6 
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10.2 
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8  points  across  vortex 


1.01 


0  5  10  15  20  25 

X 

distance  traveled:  SOOMAddlibs 

-  scale-discriminant  dissipation  preserves  structure 

-  robustness  requires  tuning  to  scheme  resolvability 

-  efficacy  limited  by  dissipation  scheme 
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(kA  y)/jt  (kA  y)/jt 


Isentropic  Vortex: 

Scale- Discriminant  Stabilization 


CFLuW~  0.01  a  =  M>  =  l 

CD04  f  1 1  Pt  B&B)  +  IF48  (2/3Pi) 


ref 

IF48  (2/3  Pi) 
EFIO 


4  points  across  vortex 


distance  traveled:  SOOMAddlibs 


-  scale-discriminant  dissipation  preserves  structure 

-  robustness  requires  tuning  to  scheme  resolvability 

-  efficacy  limited  by  dissipation  scheme 
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Conclusions/Going  Forward 


□  non-linear  instabilities  instigated  by  high  frequency  error  (i.e. 
dispersion,  aliasing) 

□  remove  error  with  respect  to  overall  scheme  resolvability  -  a 
robust  strategy 

□  scale-discriminant  dissipation  provides  stability  while  minimizing 
dissipation  error 


going  forward: 

□  non-linear  stability  (residual  filtering,  skew-symmetric  forms  etc...) 

□  incorporating  temporal  error 

□  applications  to  explicit  LES 
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